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ABSTRACT 


Until recently, there has been little literature 
available pertaining to the design of uncentralized 
squeeze film dampers (SFD). The purpose of this work 
was to design and build a rotor system supported by un- 
centralized SFD that could be operated with or without 
end seals and to conduct an experimental study of the 
system's vibrational response. Inlet pressure, viscos- 
ity and clearance were each found to have a great influ- 
ence on damper performance. Comparisons were made with 
the results found in available literature and with a 
simplified computer analysis. Suggestions are included 
for design improvements and for future work. 
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CHAPTER 1 


INTRODUCTION 

As the earth's material resources become more 
scarce and therefore more expensive, the development of 
greater sophistication in technologies which conserve 
materials becomes more economical. In the power and 
aerospace industries, more "power per pound" generally 
means lighter, higher speed, more flexible rotating 
equipment. Since all rotors have some amount of unbal- 
ance which increases with use, higher rotating speed can 
mean unacceptable vibrations and instability problems. 
The need for a damping device to control vibration and 
instability is apparent. 

Beginning with experimental use in the early 
1960's, hydrodynamic squeeze film dampers (SFD) have 
found increasing use over the past twenty years in the 
gas turbine industry. Today most modern gas turbine 
engines use SFD coupled with rolling contact bearings. 
Since rolling contact bearings have little inherent 
damping, the extra damping is necessary for safe flexi- 
ble rotor operation, operation at speeds near or in 
excess of the pi nned-pi nned critical speed of the rotor. 
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A SFD is a fluid (usually oil) filled annulus 
surrounding the bearing (or bearing housing) with a rad- 
ial clearance on the order of 10 mils. The damper bear- 
ing is "dogged" to prevent rotation but is free to 
"whirl" or process about the center of rotation of the 
rotor, while squeezing a pressure film ahead of it. The 
driving mechanism which causes the motion of the damper 
journal through the fluid film thereby creating the de- 
sired damping effect is the translational vibration of 
the rotor. 

A properly designed squeeze film damper can (1): 

a) reduce the level of forces transmitted 
through the bearings. 

b) reduce the vibrational amplitude of the rotor, 

c) smooth operation through the critical speeds. 

d) protect the rotor from sudden unbalance. 

e) protect the rotor from self-excited instab- 
i 1 i ty . 

Since rolling contact bearing life is related to the ap- 
plied forces, (a) above can mean greatly increased bear- 
ing life (a transmi ssi bi 1 i ty of 0.20 can mean an exten- 
sion of bearing life by 125 times (3)). In compressors 
and turbines, blade tip clearance is critical suggesting 
that (b) above is a clear advantage. The advantage of 
(c) is the capability of safe, flexible rotor operation. 
Catastrophic failure from blade loss during operation can 
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be avoided because of (d). Non-synch ronous , self-excited 
whirl instability can limit the rotational speed of 
rotors mounted on fluid film bearings. This instability 
can also be caused by internal friction or variable aero- 
dynamic loading with rolling contact bearings (10). In 
view of (e), SFD could be a solution to this instability 
p robl em . 

As suggested by Cookson and Kossa, SFD can be 
divided into two main categories based on basic design 
philosophy (14). The design generally used in the U.S. 
is a SFD in parallel with a flexible rotor support which 
preloads the rotor centrally in the damper bearing. The 
flexible preloading is usually accomplished with a set of 
cantilever rods referred to as a "squirrel cage". See 
Figure 1.1. On the other hand, the British design is the 
more basic uncentralized SFD with the rotor free to find 
its own position within the damper bearing. 

In the early analysis of SFD the rotor center was 
assumed to follow a circular synchronous orbit around the 
center of the damper bearing clearance circle. This ap- 
plies for centrally preloaded or vertical rotors. Though 
uncentralized SFD are simpler and less expensive, analy- 
sis of this design is non-linear and much more difficult. 
Proper design is important, however, since a poor design 
can amplify rather than attenuate transmitted bearing 
forces . (14) 
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Symmetric, Lumped Mass Rotor with 
Damped, Flexible Supports 
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A flexible support in parallel with a SFO effec- 
tively reduces the natural frequencies of the system so 
that the rigid rotor critical speeds can be traversed 
well before the operating speed. The flexible element is 
also useful for carrying the static load for a heavily 
loaded rotor. A SFD used alone does not have such a 
marked effect on the system's natural frequencies. Since 
most gas turbines now run well above the rigid rotor 
modes, the flexible element is often unnecessary (12). 

A secondary classification of SFD which applies 
to both of the above mentioned types of design involves 
the use or non-use of end seals. Non-use of damper end 
seals and low inlet pressure allows damper fluid to flow 
freely longitudinally in the damper annulus yielding a 
parabolic longitudinal pressure distribution for a short 
bearing. A damper with no end seals and a "short" bear- 
ing (L/D < .3) can be analytically approximated with the 
"short bearing solution" of the Reynolds equation. The 
fluid feed for a damper with no end seals is usually ac- 
complished with use of a ci rcumf erent i al groove in the 
center of the outer race of the damper bearing, which 
divides the damper into two parallel bearings. See Fig- 
ure 1.2. 

A damper with end seals will have no longitudinal 
flow and therefore, a constant longitudinal pressure dis- 
tribution. This type of damper can be analytically 
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approximated with the "long-bearing solution" of the Rey- 
nolds equation. Oil flow through this type of damper, 
which is required for heat dissipation, is usually main- 
tained with the use of inlet and outlet ports located 
centrally in the damper annulus 180° apart on the bottom 
and top, respecti vely . See Figure 1.3. 

Most designs include sealing devices, but the 
seals may be used to channel flow through the desired 
outlet at ambient pressure rather than to form pressuriz- 
ing end seals. See Figure 1,2, The actual sealing ar- 
rangement can take many forms, but the sealing devices 
are usually 0-rings, used in shear or compression, or 
piston rings. See Figures 1.3 and 1.4. The particular 
application typically constrains the flow rate and, 
therefore, the sealing arrangement. 

The proper design of a SFD depends on the type of 
vibration that is to be suppressed, the one type being 
the synchronous, forced, unbalance response and the other 
being nonsynchronous whirl instabilities. It is diffi- 
cult to design a SFD for all conditions in a wide operat- 
ing range. It is possible for a SFD to suppress the re- 
sponses to small normal unbalances but provide worse than 
rigid support response for large unbalances such as from 
a blade loss (10), 

When a rotor is said to be "rigid", the operating 
speed of the rotor can include the rigid rotor critical 
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speeds (cylindrical and conical “bounce" modes) but is 
kept well below the first bending critical speed (1). 
Whereas a rotor operating near or above its first bending 
critical speed is termed "flexible". 

In a linear analysis the SFD is modeled as a 
spring and damper system where the stiffness and damping 
coefficients are constant for a certain speed. In a non- 
linear analysis the damping and stiffness are functions 
of the displacement and velocity of the damper journal. 
For circular, centered, constant speed orbits the non- 
linear analysis simplifies to one of amplitude dependent 
coef f i ci ents ( 4) . 



CHAPTER 2 

A REVIEW OF SOME RECENT STUDIES CONCERNING SFO 

In 1970, Gunter published an investigation con- 
cerning the determination of the desirable values of rol- 
ling element bearing support stiffness and damping for 
rigid rotors (1). The mathematical model simulated a 
general four degree-of-f reedom unbalanced rotor mounted 
on two damped, linearly flexible supports. The stiffness 
and damping coefficients could be constant or speed de- 
pendent . 

The four simultaneous differential equations of 
motion for the rotor system were developed and then solv- 
ed with a computer for synchronous precession over a wide 
range of rotor speeds and with various support conditions. 
To simplify the analysis, isotropic stiffness and damping 
were used, yielding circular, steady-state rotor orbits. 

For optimization of support coefficients, a com- 
puter program was written that calculated the rotor amp- 
litudes and forces at each speed increment and iterated 
the speed until the maximum forces and amplitudes of rotor 
motion were obtained. For a certain stiffness, an optimum 
value of damping could be found, Gunter noted that the 
optimum value of damping was a compromise between the 

11 
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forces transmitted through the bearing supports at the 
rigid rotor critical speeds and at the maximum rotor 
speed. Additional damping attenuated forces at the cri- 
tical speeds, but raised the force level at higher 
speeds. 

A "stiff" bearing support that placed the second 
rigid body critical speed above the design speed, led to 
excessive transmitted forces with unbalance. Gunter 
also noted that any amount of damping was ineffective for 
attenuating force transmission when placed in parallel 
with stiff supports. 

Large amounts of damping were also found to be 
undesirable since excessive damping values caused force 
transmission to increase rapidly with speed. A study was 
made of the displacement phase angles of both bearings 
with respect to the radial unbalance vector. With light 
damping the phase changes at the critical speeds were 
abrupt. With moderate damping, the phase changes were 
smooth. As the rotor speed was increased, the rotor in- 
verted to rotate about its principal inertial axis. Ex- 
cessive damping suppressed the inversion of the rotor , 
thus causing high bearing forces to be generated. 

With certain support conditions there was no ob- 
servable critical speed response. Undamped, "soft" sup- 
ports could minimize rotor force transmission over the 
entire speed range considering only synchronous. 
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smal 1 -unbal ance response. However, Gunter recommended 
some damping to control instability problems (nonsyn- 
chronous whirl) and in case of shock loading, such as 
from the loss of a turbine blade. A combination of low 
stiffness and moderate damping was also found to be ef- 
fective in minimizing transmitted forces resulting from 
unbalance. With these support conditions the rotor could 
be safely operated at a critical speed. 

Gunter states that the behavior of a SFt) is 
highly non-linear so that in certain conditions a SFD can 
cause amplification of force rather than attenuation. 

In 1 973, Giberson, 4j-f- T^J-rbo Research Inc., pub- 
lished a paper acclaiming the ability of SFD to stabilize 
otherwise unstable rotors supported with journal bear- 
ings (2). When rotor speed surpasses the threshold of 
instability for the journal bearings, "whirl" amplitude 
theoretically grows without bound, possibly destroying 
the machine. Pivoted pad journal bearings can stabilize 
a system by themselves, but "pivoted pad bearings sup- 
ported by SFD is the arrangement with the greatest 
ability to stabilize a rotor bearing system that uses 
hydrodynamic bearings." 

To demonstrate the stabilizing effect of SFD, 
Giberson used a computer simulation of an analytical 
model that contained non-linear force-displacement 
(stiffness) and force-velocity (damping) relationships. 
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The equations of motion for the rotor masses were inte- 
grated through time, yielding time transient trajectories 
for the rotor discs, journal center, and damper journal 
center. 

The modeled system was a multimass rotor mounted 
on fixed pad journal bearings supported by SFD. 0-rings 
used in compression/tension acted as flexible support 
elements as well as end seals. Oil entered through an 
inlet port centered axially in the damper housing. 

Damper length, damper film cavitation pressure, oil sup- 
ply pressure, and rotor speed were kept fixed for all 
presented simulation runs. The objective was to "tune" 
the system by varying three support parameters; the rad- 
ial damper clearance,. 0-ring stiffness, and oil viscos- 
ity. Decreasing the clearance, raising the viscosity, or 
stiffening the 0-rings tended to increase support rigid- 
i ty . 

The results showed that the system was unstable 
with rigid supports and was stable for only one out of 
nine choices of parameters. As Giberson points out, SFD 
are capable of solving whirl problems, but they must be 
"tuned" (tested to determine the correct combination of 
design parameters for the particular system) to be effec- 
tive. In fact, an improper design can lower rather than 
increase the threshold of instability. 
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In 1974 Mohan and Hahn published a paper on the 
"Design of SFO Supports for Rigid Rotors" (3), The pub- 
lication includes a general design guide for centrally 
preloaded SFD, supporting a rigid rotor, mounted on rol- 
ling contact bearings. The authors also investigated the 
transient solutions of the fully non-linear equations of 
motion to discover the effect of central preloading. 

The rotor that was modeled was a rigid, symmetric 
rotor mounted on ball bearings, supported by retainers 
springs and SFD with no end seals. The fluid film forces 
necessary for solving the equations of motion for the 
system were obtained with the use of the short bearing 
approximation of the Reynolds lubrication equation for 
constant fluid properties. The assumptions used in solv- 
ing for the film forces required that; 

1) an incompressible lubricant and an average 
fluid viscosity in the squeeze film be used, 
(constant properties) 

2) fluid inertia and turbulence effects be neg- 
ligible. (necessary to derive Reynolds equa- 
tion) 

3) the length to diameter ratio (L/D) be less 
than 1/4. (short bearing) 

4) the end pressures be ambient, (no end seals, 
low feed pressure) 
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5) the tensile forces in the film be negligible 
despite below ambient pressures. 

6) the extent of the positive film pressure re- 
gion be the same for a full film and a cavi- 
tated film. 

The Reynolds equation was integrated over the 
length of the bearing to yield an expression for film 
pressures in terms of damper journal position and velo- 
city. The expression for the pressure was then integrat- 
ed over the positive pressure region of the damper, 
yielding the desired tangential and radial film force 
components. The equations of motion for the rotor system 
could then be solved for each time step, yielding the 
journal position and velocity and the fluid film forces. 

The solutions were dependent on the following 
non-dimensional parameters: 


1) 

B = 

uRL3 

(beari ng parameter) 

2) 

U = 

u 

msC" 

(unbalance eccentricity) 

3) 

0) 

IT 


(ratio of operating speed to 
design speed) 

4) 

Ws 

(jJ 2 

R“ 

(spring parameter) 
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5) (gravity parameter) 

6) spring preload parameters 

Consideration of the transient, non-linear equa- 
tions of motion yielded one, two, or three stable modes 
of operation (at design speed above rigid critical speed) 
depending on the operating conditions, 

1) A circular synchronous journal orbit called 
the "inverted mode" (the rotor rotates about 
its principal inertial axis, the journal dis- 
placement and unbalance force are out of 
phase by 180°). This is a desirable mode in 
that it is characterized by a small journal 
orbit and low transmi ss ibi 1 i ty (T). 

2) The "whirl mode" characterized by a circular, 
synchronous orbit but with an uninverted 
rotor (the journal displacement and unbalance 
force are out of phase by 110° or less). 

This mode is undesirable in that it has a 
large journal orbit (high eccentri ci ty ) and a 
T > 1. 

3) An extremely undesirable mode with a nonsyn- 
chronous, noncircular orbit including violent 
excursions, high eccentricity (e) and high 

t ransmi ss i bi j i ty , 
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The third mode was only observed for low values 
of the bearing coefficient (B < 10"2) with certain bound- 
ary conditions. Operation in either the first or second 
mode depended on the values of the damping and unbalance 
parameters (B,U), Note the agreement with Gunter's 1970 
analysis concerning the inverted and uninverted modes. 

While central preloading eased the steady state 
solution by reducing the number of parameters involved, 
it did not adversely affect the transmi ss i bi 1 i ty or oper- 
ation modes. The design data found in this paper is re- 
stricted to steady state conditions for centrally pre- 
loaded SFD. Linear theory was found to apply adequately 
for low unbalance (U < .2) for obtaining the frequency 
response from knowledge of the resonant speed and trans- 
missibility. For U > .2 the system had to be treated as 
non -1 i near. 

When U was plotted against steady state T for a 
certain bearing it was found that for increasing U there 
was a critical unbalance (U;-) at which T suddenly in- 
creased (jumped). The "jump phenomenon" occurred as the 
system changed from operation in the inverted mode to 
operation in the uninverted mode. For high unbalances, 
the rotor did not invert at all. 

It was found that lowering the spring rate had 
the desirable effects of raising Uc and lowering T at 
both the design speed and resonant speed. This suggests. 
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as Gunter did in 1970, the use of as low a spring rate as 
possible for maintaining central preloading. 

It was also discovered and shown graphically that 
a large B is desirable at resonant speed but a small 
value of B is desirable at design speed (above critical) 
from the standpoint of producing a low T, This suggests 
(as Gunter's analysis did in 1970) a compromise value of 
B that allows for safe acceleration through resonance but 
low T at design speed. The authors suggest that for an 
uncertain unbalance condition, it is best to choose a 
value of B that maximizes the critical unbalance. Plots 
are provided for designing SFD with flexible supports for 
centrally preloaded rigid rotors. 

In 1975 Cunningham, Fleming, and Gunter published 
an analytical investigation on the design of a centrally 
preloaded SFD for a multimass flexible rotor mounted on 
rolling contact bearings (4), The paper demonstrates a 
technique for the use of single mass, symmetric, flexible 
rotor analysis to optimize the stiffness and damping for 
a symmetric five mass rotor. The single mass analysis 
was taken from a 1972 NASA Report by Kirk and Gunter. 

To aid in the analysis of the single, lumped mass 
flexible rotor system, the following dimensionless para- 
meters were considered: 


_ bearing housing mass 
rotor mass 


(mass ratio) 
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K 


support stiffness 
rotor-bearing stiffness 


(stiffness ratio) 


total support damping 
rotor-bearing damping 


(damping ratio) 


The flexible supports effectively reduce the 
first critical speed to some fraction of the pinned-pinn- 
ed rigid support, bending critical speed. Increased sup- 
port damping was found to decrease the rotor amplitude at 
the critical speed but increase the amplitude at the 
rigid support critical speed. In other words, increased 
damping tended to "lock up" the rotor support. The "op- 
timum" value of damping provided about equal rotor ampli- 
tude at the critical and rigid support critical speeds. 
(Note the increased complication due to the flexible 
rotor). It was demonstrated that for high rotor speed 
(above rigid support critical), T increased with increas- 
ing damping while rotor orbit amplitude was unaffected, 
as noted for rigid rotors in the previous paper. In 
choosing the damping, the forces and amplitudes of motion 
throughout the operating range must be taken into consid- 
eration. 

For a given mass ratio (M), rotor amplitude of 
motion is minimized for a tuned system (K = M). The 
smallest possible value of M is desirable for lowering 
rotor amplitudes of motion. However, too small of a 
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value of K leaves the system susceptible to shock and 
self-excited whirl problems (in the absence of proper 
dampi ng ) . 

For a multimass flexible rotor, if the maximum 
rotor speed is below the second critical speed, single 
mass theory is shown to be adequate in some cases. The 
rotor that was simulated was a five mass symmetric rotor 
mounted on ball bearings with damped, flexible supports 
that were to be designed for optimum performance. Oil 
feed to the SFD was through a central ci rcumferenti al 
groove. Piston rings and metered orifices controlled oil 
flow. The rigid support critical speeds were first cal- 
culated with computer simulation. For support design 
calculations, the five masses were replaced by a single, 
centered mass that resulted in the same first critical 
speed. The optimum support damping was found for this 
single mass rotor configuration. See Figure 1.1. 

The assumption of a centralized circular orbit of 
the damper journal was used in designing a SFD that would 
provide the desired damping. Dimensionless stiffness and 
damping coefficients were obtained from the short-bearing 
approximation of the Reynolds equation assuming full 180® 
cavitation, no tensile stress and constant properties, 

r , - - 2 = 
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Note the non-linearity in terms of orbit size (eccentri- 
city). Damping can be seen to be a function of damper 
clearance, radius, length, oil viscosity, and damper 
journal orbit. 

A maximum eccentricity of e = .4 at the critical 
speed was chosen for the design, somewhat arbitrarily. 

The maximum support amplitude of motion should be about 
equal to the mass eccentricity (rotor mass C.G. offset 
from axis of rotation) for optimum damper performance 
according to the authors. Thus, from knowledge of the 
maximum expected unbalance, the damper clearance could he 
calculated. Considering the radius and viscosity to be 
given, the desired damping then determined the damper 
1 ength . 

An unbalance response computer program was used 
to discover the effectiveness of the designed supports 
for the five mass rotor. The dampers successfully atten- 
uated rotor orbit amplitudes of motion and bearing forces 
at the first critical speed by factors of 16 and 36, re- 
spectively, from the rigid support case. The vibration 
and forces were also reduced at higher critical speeds. 
However, for unbalances much greater than the design un- 
balance, the SFD supports provided worse performance than 
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rigid supports. 

In 1979 Tonneson, of the Technical University of 
Denmark, published a study on experimental squeeze film 
bearing orbits (5). The purpose of the investigation was 
to compare experimentally measured damping coefficients 
with those obtained from simple linearized theory. Comp- 
arisons were made for both eccentric (offset) and concen- 
tric (centralized) damper journal orbits. (See Figure 
2.1 ) . 

Two typical experimental apparatus are described. 
The first was a ball bearing supported rigid rotor mount- 
ed on SFD. The oil feed was accomplished with a circum- 
ferential groove in the damper land. End flow seals were 
provided. The length to diameter ratio was 0.34. The 
second test rig used was a similar SFD but the bearing 
housing was excited by an electrodynamic shaker. 

The first apparatus used a flexible support, with 
a differential adjustment for centering the damper jour- 
nal, in parallel with the SFD. Both rigs were instrumen- 
ted for measurement of forces, displacements, phase 
angles, frequencies, and extent of cavitation. 

The theory used for comparison with experiments 
was the short-bearing approximation of the Reynolds equa- 
tion for small, circular, centered orbits (e < .3 and 
constant). The damping coefficients for this linear (in- 
dependent of position, e, and velocity) approximation are 
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Bd = Bj(-By=^yR 


for a 180° positive pressure film (expression should be 
doubled for a 360°, uncavitated film). 

The results for the concentric orbit case show 
that the measured damping coefficients were between the 
theorized damping coefficients for 180° and 360° film. 

For the first test rig, for a constant viscosity and 
clearance, the damping coefficients were shown to be in- 
dependent of frequency until the orbit amplitude exceeded 
a certain value. The shaker driven test rig yielded data 
showing that varying the viscosity, radial clearance, or 
orbit size yielded a linear relationship of with the 
frequency for small amplitudes. 

The experiments also revealed that the end seals 
added stiffness and damping which could be significant 
for low values of 8,^ for the SFD. The theory assumed at- 
mospheric end conditions (short bearing), but the test 
rigs used end seals which change the end conditions to 
that of the inlet pressure so that the long bearing ap- 
proximation may have been more valid. The supply pres- 
sure was found to affect the film extent and therefore 
the damping. 

The author suggests that the eccentric case 
(damper orbits about a position off-center of the damper 
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housing) is more typical of SFD operation since exact 
central preloading is difficult due to built in toler- 
ances. Reasonable agreement with theory was again found 
for small offsets and small orbit amplitudes. For large 
offsets the damper characteristics were pronouncedly non- 
linear, for instance, lowering the viscosity at a large 
offset drastically increased the orbit amplitude. 

For design applications, Tonneson suggests that a 
small L/D (about .1 to .2) is desirable and C/R should he 
no less than .002. A soft spring support is again recom- 
mended for effective SFD performance. For linear theory 
to apply, the flexible support must be stiff enough to 
accomplish approximate central preloading. 

In 1977 Bansal and Hibner of Pratt and Whitney 
published an "Experimental and Analytical Investigation 
of (dynamic) SFD Forces Induced by Offset Circular Whirl 
Orbits." (6) In agreement with Tonneson, the authors 
suggest that offset orbits are more realistic for flex- 
ible damped supports. The offset is the distance between 
the steady state orbit center and the damper housing 
center, not to be mistaken with the eccentricity, e, which 
is the instantaneous distance between the journal center 
and housing center (see Figure 2.1). Offset circular 
orbits were induced mechanically on a single SFD support- 
ing a rolling contact bearing. The experimentally measur- 
ed pressure distribution was used to calculate the dynamic 
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forces. The purpose of the investigation was to study 
the effect of varying the inlet pressure, viscosity, 
journal speed and damper eccentricity on damper forces. 

A motor driven eccentric cam produced the requir- 
ed offset circular orbits (actually, the dynamic orbits 
were elliptical because of the dynamic forces of the 
damper and the tolerances in the roller bearing). The 
damper oil inlet was a pressure fed ci rcumf erenti al 
groove at one end of the damper annulus, which provided 
uniform flow into the damper, and the oil exit was a 
single port on the other end of the annulus. 0-ring end 
seals were used. Eight pressure transducers, axially 
centered and distributed circumferentially about the 
damper housing, sensed the film pressure and proximotor 
probes sensed the orbit shape and size. 

The analytical solution for the pressure profiles 
was produced with the use of the long bearing approxima- 
tion of the Reynolds equation assuming constant proper- 
ties. The experimentally measured values of orbit posi- 
tion and shape were used in the analytical solution and 
experimental data reduction. The dynamic pressure at any 
circumferential angle, 0, P(0), within the damper annulus 
was determined by superimposing the effect of the inlet 
pressure , P-jn» op The hydrodynamic pressure Ph(0); 

P(9) = Ph( 9) + Pin* P(9) for 0 from 8^ to was 
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less than the cavitation pressure, of the oil, then 
the film had cavitated. Note that Pin could be increased 
so that P(9) > pQ for all 0 yielding an uncavitated film. 
The cavitated portion could range from 180® (Pin * 
bient pressure) to 0® (full film). When cavitation was 
present, the damper film stiffness and damping were very 
sensitive to a change in the inlet pressure. For circu- 
lar centralized orbits, an increase in inlet pressure in- 
creased the damping forces and decreased the stiffness. 

For a centered circular orbit, the eccentricity 
is constant and the circumferential pressure distribution 
measured at any point of the annulus is identical, so, 
dynamic pressure measurement requires only one trans- 
ducer. For an offset orbit, the dynamic pressured induc- 
ed by the orbit is a function of the instantaneous posi- 
tion of the damper center and many pressure sensors are 
necessary. The pressure measured by each transducer 
would be dependent upon damper position and the location 
of the transducer on the rig. 

The data reduction yielded an instantaneous ex- 
perimental pressure profile for every 30® position of the 
damper center. Each profile was integrated to find the 
radial and tangential force for each damper position. 

The predicted pressure distributions, including inlet 
pressure effects and measured orbit parameters, were 
within 12% of the measured pressure distributions. The 
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comparisons of forces showed much poorer correlations. 

The authors suggested that small errors in measuring or- 
bit parameters could cause significant errors in the ana- 
lytical results due to the complex data reduction. 

As noted previously the effect of inlet pressure 
is important with a cavitated film. The general trend 
was noted to be a decrease in radial force and an in- 
crease in tangential force for an increase in inlet pres- 
sure which the authors note is the same trend as for 
centralized orbits (increase in P^n means an increase in 
damping and a decrease in the stiffness of the film). 

In 1977 Hibner, Bansal and Ruono of Pratt and 
Whitney published a paper on the control of instability 
for intershaft SFO (7). An earlier analytical and ex- 
perimental investigation had shown the existence of an 
intershaft viscous damper instability. The intershaft 
location in question was between the low and high speed 
shafts of a gas turbine engineer, at a roller bearing 
location. The outer race of the roller bearing was fixed 
to the outer, higher speed shaft, and the inner race was 
fixed to rotate with the inner, lower speed, fan shaft. 
The damper annulus was located between the inner race 
housing and the low speed rotor, so, the SFO journal and 
housing rotated with the low speed rotor. It was found 
that the "intershaft SFO would excite an instability at 
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low rotor speed as the cantilever natural frequency of 
the high speed rotor was approached.'* 

The authors demonstrated that a flexible element 
in parallel with the SFO could add the necessary radial 
stiffness needed to raise the instability threshold above 
the operating speed. 

In 1979 Hahn published an analytical investiga- 
tion on the "Unbalance Behaviour of Squeeze Film Support- 
ed Rigid Rotors." (9) Equilibrium load capacity and 
t ransmi ssi bi 1 i ty data for a wide range of operating con- 
ditions are presented for centrally preloaded and unload- 
ed vertical rigid rotors mounted on rolling element bear- 
i ngs , 

Hahn noted that the performance of a SFD is 
strongly influenced by: 

1) the damper design parameters 

2) the degree of rotor unbalance 

3) the damping inherent in the system 

4) the operating speed 

5) the degree of oil supply pressurization. 

The theoretical model of the damper came from the 
short bearing approximation of the constant property Rey- 
nolds equation. The analysis was restricted to steady 
state, centralized, circular, synchronous orbits. The 
solution for the orbit eccentricity was presented 
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as plots of two speed' dependent dimensionless parameters 
A vs. B/U. 

A = (1 - l/fi2)/u 

B = (Load Number) 

6 = — 

«c 

U = Unbalance Parameter 

Two plots are given, for an unpressurized and a fully 
pressurized film. Given a rotor system (lOc.m.U), the SFD 
design parameters (y,R,C,L), and the operating speed it 
is possible to locate predicted values of eccentricity 
and the system's operating line on the plots. The op- 
erating line would indicate the speed at which the maxi- 
mum vibration amplitude would occur. 

On the orbit eccentricity plots for unpressurized 
bearings there is a region where a point on the plot can 
suggest three possible equilibrium orbits (note the 
agreement with the earlier papers of Hahn and Gunter). 

The intermediate orbit has been shown to always be un- 
stable, leading the author to label this region the bi- 
stable region. During operation in the bistable region 
it is possible for the system to jump from one to the 
other stable orbit. The possibility of jump depends on 
the conditions at entry into or exit from the region. To 
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avoid a jump to a higher orbit (higher eccentricity is 
associated with higher transmi ssi hi 1 i ty ) , it is. best to 

maintain operating conditions that keep the system out of 
the bistable region. 

Unpressurized Dampers 

Hahn recommends 6 > 2 for maintaining T < 1 and 
operation generally above the bistable region on the de- 
sign plots. There is significant vibration attentuation 
possible for U < .43. For unloaded vertical rotors, there 
is no advantage in using external flexible supports in 
parallel with the SFD. 

Fully Pressurized Dampers 

There is no bistable region, but external radial 
support stiffness in parallel with the SFD is required 
for the stability of the solution. In this case, verti- 
cal rotors must have retainer springs. For 6 > 2, T 
will always be less than unity despite the value of un- 
balance, but as U gets large, T approaches one. 

Note that vibration isolation can only be achiev- 
ed at the expense of increasing rigid rotor vibration, so 
rotor motion constraints (tip clearances! must be consid- 
e red. 

In 1979 Vance of Gas Turbine Laboratories pub- 
lished a review of rotor dynamics which included the cur- 
rent role of SFD as a solution to dynamical problems (10) 
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The author stated that it was still not clear that SFD 
could be designed purely analytically for safe reliable 
performance for all cases. 

Vance pointed out that the then current analyti- 
cal techniques used the constant property Reynolds equa- 
tion which does not apply for certain practical cases. 

The questionable assumptions generally used are; 

1) negligible fluid mass inertia effects 

2) single phase fluid film (no gas bubbles in 
the pressure region), 

3) constant vi scos i ty throughout the film. 

The author also pointed out that the Reynolds 

equation is highly nonlinear (despite the typical use of 
a linear analysis) and that designers must be aware of 
instability problems, the neglect of which has been a 
costly error in the past. 

U.S. army turboshaft helicopters now run at 
" supercri ti cal " speeds (above the third critical flexural 
speed) using SFD to attenuate vibration and suppress 
whirl instability. Journal bearings are not used in air- 
craft engines because of potential catastrophic failure 
if oil flow is blocked. Vance argued that SFD should not 
be used to suppress instabilities in aircraft engines for 
the same reason. 

In 1979 Hibner and Bansal published an experimen- 
tal and theoretical investigation into "The Effects of 
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Fluid Compressibility on Viscous Dampers," (11) The in- 
vestigation was begun because of the questionable valid- 
ity of the Reynolds equation when cavitation exists in 
the squeeze film. The effects of cavitation could ex- 
plain why correlations between theory and experiment for 
SFD have been excellent in some cases and poor in others. 

The experimental test rig consisted of a control- 
led orbit, end sealed SFD, Circular, centered damper 
journal orbits were produced with an offset cam. A 
single pressure transducer measured the circumferential 
pressure distribution and probes measured the orbit 
eccentricity. 

The theoretical analysis made use of the conven- 
tional long-bearing approximation of the constant pro- 
perty Reynolds equation. The experimentally measured 
dynamic values of eccentricity were used in the theoreti- 
cal calculations. 

The results showed better theoretical/experimen- 
tal correlation for a decreasing extent of cavitation. 

For a cavitated film, theory overpredi cted the pressure 
distribution. Raising the supply pressure or lowering 
the speed decreased the extent of cavitation and improved 
the correlation. With a cavitated damper, the supply 
pressure could not be superimposed on the hydrodynamic 
pressure, since an increase in supply pressure yielded an 
increase in hydrodynamic pressure. 
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It was hypothesized that the "deviation from 
theory is due to fluid compressibility caused by sub-at- 
mospheric cavitation pressures liberating dissolved gas 
and creating a two phase fluid." The liberated bubbles 
"persist within the pressure region thereby reducing the 
hydrodynamic pressure." An increase in the supply pres- 
sure yielded less cavitation, less bubbles, and higher, 
closer to predicted hydrodynamic pressures. 

The effect of compressibility when ignored in a 
design can yield improvement or degradation in the de- 
sired vibrational response of a SFD. The authors suggest 
designing SFO in the range where the incompressible Rey- 
nolds equation works well or do experimental testing of 
the design because the influence of a two-phase fluid on 
a viscous damper is complex. 

In 1979 Holmes of the University of Sussex, U.K., 
published a study on the control of rotor vibration using 
SFD with and without a parallel flexible element (12). 

The experimental and analytical investigation includes a 
general design philosophy. 

SFD with a Flexible Support 

A SFD as used in a typical application provides 
light, non-linear damping which makes it susceptible to 
two classical disadvantages of such a damper, jump pheno- 
menon and subharmoni c resonance instability. Holmes 
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notes that, as earlier shown by Hahn, jump can only occur 
when cavitation is present, but cavitation is not nec- 
essary for subharmonic resonance to occur. 

The experimental test rigs had rotors (one rigid, 
one flexible) mounted on journal bearings supported with 
SFD and centrally located with adjustable springs. 

The results from the rigid rotor test rig showed 
good agreement with conventional Reynolds theory for a 
fully pressurized, uncavitated damper. Note the agree- 
ment with Hibner and Bansal on this point. When cavita- 
tion is present, the author suggests that the uncertainty 
of the extent of the cavitation makes theoretical predic- 
tions difficut. Holmes recommends that the supply pres- 
sure should be as high as possible for the particular 
design to reduce the extent of cavitation. 

For the flexible shaft test rig the optimum 
squeeze film damping for reducing the maximum shaft amp- 
litude at the first and second critical speeds was found 
to be about the same for both theory and experiment. 

For light squeeze film damping (below the optimum) the 
system was stable and well damped in the entire operating 
speed range, but for large amounts of damping instability 
thresholds were observed. For moderate damping (near the 
optimum value) steady sub-synchronous whirling would ap- 
pear at a certain speed, grow and die out as speed was 
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increased. The whirl was indicative of the non-linearity 
of the damper and journal bearing oil films. 

Linear analysis can be used for SFD with parallel 
flexible supports for damper amp litudes less than one 
third of the radial clearance. Above this. Holmes sug- 
gests a quasi-linear analysis using amplitude dependent 
stiffness and damping coefficients, 

SFD Used Alone as Vibration Isolators 

The lift for this type of damper emanates from 
the non-linear effects of the oil film. Jump phenomenon 
and subharmonic resonance are also possible for this type 
of support. Generally, theoretical predictions must be 
performed with non-linear analysis. For regions near the 
critical speeds where the dynamic load is much greater 
than the static load, circular, concentric orbits can be 
assumed and quasi-linear analysis performed. The extent 
of cavitation greatly affects the static and dynamic load 
carrying capacity. In fact, the assumption of an uncavi- 
tated film leads to an elimination of the load carrying 
capacity. This appears to be in agreement with Hahn's 
statement that flexible supports in parallel with the SFD 
are necessary for stability for pressurized dampers. 

The experimental test rig consisted of a rigid 
rotor mounted on rolling contact bearings supported by 
SFD with no end seals. A negative pressure region (about 
-40 psi ) was measured suggesting that tension forces may 
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exist. The measured pressure peaks and dynamic forces 
were less than predicted by a conventional short-bearing 
solution of the Reynolds equation. Holmes, in agreement 
with Bansal and Hibner's previous paper, suggests that 
cavitation bubbles trapped in the oil film are respon- 
sible for the lower pressure peaks. Since the journal of 
a SFD does not rotate, there is no inducement for the 
bubbles to leave except the supply pressure. The author 
showed that by varying the cavitation extent, the theore- 
tical model for the pressure distribution could be ad- 
justed to agree with experimentally measured pressures. 

In 1980 and 1981, Cookson and Kossa published the 
results of their experimental and analytical investiga- 
tions into the effectiveness of SFD, without parallel 
flexible supports, used for supporting flexible rotors 
(13,14). 

The short-bearing approximation of the Reynolds 
equation was used with a 180° film model, thus assuming 
atmospheric inlet and outlet oil flow conditions (no end 
seals, low inlet pressure). In developing the equations 
of motion, the rotor was assumed to be symmetric, with 
its mass lumped at the bearing stations and at the middle 
of the rotor. The theoretical analysis required constant 
speed and that the rolling contact bearings introduced no 
exciting forces (no source for generation of whirl in- 
stability). 
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Four simultaneous non-linear differential equa- 
tions were developed as the equations of motion for the 
system with the damper journal and rotor disc displace- 
ments and velocities as the independent variables. The 
solutions of the equations were found to be dependent on 
six independent dimensionless parameters three for the 
damper and three for the flexible rotor system. The 
equations were integrated forward in time for certain 
values of the system parameters yielding steady state 
solutions. 

The experimental test rig consisted of a disc 
centrally located on a rotor mounted in rolling contact 
bearings supported by SFD. Interchangeable damper parts 
produced a wide range of bearing and gravity parameters; 



A circumferential groove in the damper housing provided 
the oil inlet. Oil supply pressure was kept low (about 2 
psi) and no end seals were used. Noncontacting trans- 
ducers measured the disc and journal orbits. 

The predicted orbits and the measured orbits were 
found to be similar. Disc amplitude vs. frequency plots 
showed a smooth transition through the critical 
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speed, not a "peaky" response. The maximum amplitude 
still occurred approximately at the pi nned-pi nned criti- 
cal speed. The authors suggest that the dampers operate 
by allowing for bearing movement instead of disc movement 
(for a flexible rotor). Remember that, as Hahn noted, 
for a rigid rotor the price of vibration isolation is 
shaft movement. 

The theory was reasonably accurate in predicting 
the experimental results for the response curves and for 
the orbits for various parameter combinations. The 
authors suggest that a better film model would have im- 
proved the correlation, 

A vertical rotor could be analyzed by setting the 
gravity parameter equal to zero. For this case, the 
rotor orbits become circular and concentric as noted by 
Holmes for rotors with light static loading. 

The orbit amplitudes were found to vary little 
for a wide range of values of the bearing parameter, sug- 
gesting that the land length is not an important para- 
meter. In general, for greatest damper efficiency, 
Cookson and Kossa recommend that; 

1) The mass ratio, M. should be kept small. 

This allows for freer movement of the damper 
journals. Note the agreement with the 1975 
paper by Cunningham, Fleming, and Gunter for 
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a SFD with flexible supports supporting a 
flexible rotor, 

2) the bearing parameter should be approximately 
equal to ,1, 

3) the gravity parameter should be less then ,1 
(light static/dynamic load ratio). This also 
allows for freer movement of the damper jour- 
nals. 


The authors recommend a transient analysis for shock 
loading situations. 




CHAPTER 3 


EXPERIMENTAL SET-UP 

The rotor system was designed for use on a 
Bently-Nevada rotor dynamics test rig. The rotor was 
powered by a 1/10 horsepower infinitely controllable 
drive motor capable of driving the test rotor at speeds 
in excess of 10,000 RPM. The system included a shaft 
supported by two sets of preloaded duplex ball bearings, 
with each duplex set mounted in an uncentralized SFD. A 
quill shaft was used to couple the motor and rotor. 

Rotor discs of various weights could be located at any 
position between the bearing stations. See Figure 3.1 

The preloaded duplex ball bearings provided for 
stiff bea ri ng-rotor coupling, so that the bearings could 
be assumed to be rigid in comparison with the rotor and 
damper film. The preload on the bearings also increased 
the rolling friction which decreased the maximum attain- 
able rotor speed and increased heat generation. 

The ball bearing housings which preloaded the 
duplex bearings acted as the damper journals. The damper 
journal diameters were 3.3125 inches; the journal land 
lengths were .475 inches. So, the length to diameter 
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Figure 3.1 


Experimental Test Facility Rotor 
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ratio was .143, suitably short for a short bearing ap- 
proximation of the Reynolds equation to apply. 

The damper housings had interchangeable inserts 
allowing the radial clearance for the fluid annulus be- 
tween the damper bearing and journal to be varied from 
.004 to .020 inches. Sealing between the jou rnal and 
damper bearing was accomplished with 0-rings in shear. 

The lateral clearance between the damper journal and 
housing was dimensioned to yield .007 inches of "squeeze" 
on each 0-ring. This amount of 0-ring squeeze was deter- 
mined with a trial and error approach to he near the op- 
timum, sealing adequately under the required oil pressure 
yet contributing minimal impedance to journal motion. 

C i rcumf erenti al grooves with outlet ports were located in 
the sides of the damper housings at the edges of the 
fluid film annulus. Ports were also provided to check 
the static oil pressure in the grooves at the end of the 
annulus with a pressure gauge. Also, inlet and outlet 
ports were located centrally in the damper inserts and 
housings at the top and bottom of the dampers. Flow 
metering valves at the outlet ports allowed for the end 
seal or no end seal conditions (or any amount of end 
leakage or end pressure desired). On/off valves at the 
inlet ports allowed for circumferential groove or port 
oil feed into the damper annulus. With the appropriate 
ports closed or open the damper configuration could he 
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modeled validly by either the short or long bearing ap- 
proximation of the Reynolds equation. An antirotation 
pin was threaded through the outside of the damper hous- 
ing, through the inserts and into an oversized hole in 
the journal, allowing journal translation but not rota- 
tion. See Figure 3.2. Also see Appendix B for drawings 
of damper parts. 

The steel rotor shaft was 13/16 inches in dia- 
meter and 12 inches long between bearings. The quill 
shaft used as the flexible coupling between the motor and 
the rotor was 6 inches long and an 1/8 inch in diameter. 
See Appendix A for rotor and quill sizing calculations. 
The flexible coupling eliminated the need for exact 
alignment of the motor and the two rotor bearings. Two 
rotor discs could be located at any position between the 
bearings with set screws. Each rotor disc had threaded 
holes through its thickness located symmetrically about 
its outside circumference for the positioning of balanc- 
ing weights. The discs were modeled after those that 
were furnished with the Bently-Nevada test rig. With 
both discs axially centered on the shaft, the first flex- 
ible mode could be excited. 

A gear pump was used to pump the oil from a heat- 
ed bath through the SFD. A metered bypass line provided 
a wide range of damper inlet pressure and flow rate. 
Heating the SAE 1 OW oil provided a wide range of damper 
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Figure 3.2 Schematic Side View of Experimental Damper 
Bearing Configuration 
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fluid viscosity. Copper-Constantan thermocouples were 
used to monitor the oil inlet and outlet temperatures 
from the dampers. A thermalwell was used for monitoring 
the pressurized inlet temperature. 

Vibration detection was accomplished with Bently 
Nevada non -contact i ng proximity displacement transducers 
The output of the proximitors was wired into a digital 
vector filter (DVF2) also manufactured by Bently-Nevada. 
The DVF2 provided a digital readout of peak to peak dis- 
placement, RPM, and phase angle for location of the 
rotor's "high spot" for balancing. Two Tektronix oscil- 
loscopes, an X-Y-Y plotter, and an HP spectrum analyzer 
were also useful aids in studying and recording vibra- 
tional response. The spectrum analyzer was especially 
useful for investigating nonsynch ronous responses. A 
five thermocouple input, digital readout potentiometer 
was used to monitor oil temperature into and out of the 
SFD. See Appendix B for the wiring diagram. 



CHAPTER 4 


EXPERIMENTAL METHOD AND PURPOSE 

The initial purpose of the experimentation was to 
verify that the test apparatus could work properly as 
designed and that operational data could be taken as 
desired. Once the system was debugged the next step was 
to conduct a brief parametric study. 

The parameters that could be varied for both the 
"short" and "long" bearing configurations of the damper 
were as foil ows : 

1 ) rotor speed 

2) unbalance 

3) clearance 

4) oil temperature 

5) oil inlet pressure 

The second purpose was then to discover the effect on the 
rotor system's performance of varying each parameter 
while holding the others constant. Comparisons could 
then be made with other authors results and with a sim- 
plified analysis and conclusions made. 

The third and final purpose was to recommend pos- 
sible improvements to the test equipment and to suggest 
further work. 
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Bode plots of peak to peak amplitude and phase vs. 
rotor speed could be made at both damper locations and at 
any point along the shaft for both x and y positions. 

Mso, damper and shaft orbits could be obtained and 
photographed for a particular rotor speed. The orbits 
(oscilloscope tracings) have a gap in them made by a 
reference mark useful for phase studies. The phase 
mentioned here is the angle between the high spot of the 
rotor and a reference mark on the shaft. The meaningful 
information is the change in this angle with speed (de- 
noting changes of angles between the high spot and the 
unbalance position). 

The first task was to balance the rotor after the 
system was assembled with the nul 1 -cl earance inserts in the 
dampers providing "rigid" bearing supports. The nVF2 aided 
in this endeavor. Figure 4.1 shows the midspan vibrational 
response of the balanced rotor. Note that all Bode plots 
were made from a location two inches from the midspan loca- 
tion. Comparing with the other plot on the same figure for 
a small unbalance, it can be seen that the residual unbal- 
ance is small. The sharp critical speed response on Figure 
4.1 and the abrupt phase change shown on Figure 4.2 for a 
small unbalance demonstrates the high stiffness and lack of 
damping of the ball bearings. The first critical speed of 
the rotor can be seen to be about 9200 RPM. 
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Figure 4.1 



Figure 4.2 
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In practice, the four damper mounting screws on 
each SFD were used to adjust the 0-ring squeeze to seal 
adequately, yet contribute minimal impedence to journal 
motion. Shims were then required to restrain the damper 
bearings from vibrating relative to their respective 
housi ngs . 

When damper orbit studies were made the GAP 
button on the DVF2 was used to obtain the approximate 
orbit center coordinates within the clearance circle. 

The GAP button provides a measure of the distance between 
the proximotor probes and the rotor. This method was 
also used to measure the radius of each clearance circle. 

Figure 4,3a and b are examples of typical photo- 
graphs of the damper and shaft orbits, respectively, 
taken from oscilliscope tracings. The rotation of the 
shaft is clockwise with respect to all orbits presented 
in tnis paper. Copies of each orbit photograph of inter- 
est were reduced to their proper scale, and placed in the 
correct position and orientation on drawings of their 
particular clearance circle. This was done because know- 
ledge of the location (offset from center) of the orbit 
within the clearance circle was found to be useful. 
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b. Shaft Orbit, same as above except 2 psi, 2.5 MIL/DIV. 


a. Damber Orbit, Groove Inlet, c - .01", 80°F, 9000 RPfl, 
10 psi, 3.0 gm unbalance, 1 MIL/DIV. 


Figure 4-3. 


Photographs of Oscilliscope Tracings of 
Damper and Shaft Orbits 


CHAPTER 5 


RESULTS 

The experimental results are presented here for 
the rotor system supported by SFO without end seals (an- 
nular oil i nl et/outl et ) and with end seals (port inlet/ 
outlet). 

SFD Without End Seals 

Figure 5.1, a superpositioning of four separate 
damper orbits, shows the effects of changing the unbal- 
ance while holding the other parameters constant. In- 
creasing the unbalance increased the vibration of the 
rotor which resulted in a further "lifting off" of the 
dampers from the bottom of the clearance circle as well 
as larger elliptical orbits. 

Figures 5.2 through 5.6 show amplitude plots of 
the rotor's x-axis for different values of unbalance. 

Note that the trends are consistent despite the clearance 
size. All cases show that at least as a small amount of 
effective damping was present, yielding a response super- 
ior to the "rigid" response of Figure A.l. This is appar- 
ent from the lower amplitudes and the less abrupt phase 
shifting with additional dissipation. As more unbalance 
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X 1000 RPM 


Figure 5.2 
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Midshaft Phase vs. Rotor Speed 

Groove Feed, 80°F, C = .010", = 2 psi 

1 U = .0662 (1 gm) 

2 U = .1987 (3 gm) 



Figure 5.3 



Figure 5.4 






Pn'ASE 


Midshaft Phase vs Rotor Speed 


Groove Feed, 80®F, C = .006", P.^ = 2 psi 

1 u = 1 gm 

2 u = 2 gm 



Figure 5.5 
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was applied, larger orbits were observed and the effec- 
tive damping increased, displaying the amplitude depen- 
dent damping characteristics of the SFD. As noted by 
Cookson and Kossa for this type of dampers, the critical 
speed did not change appreciably. As will be described 
later, damping tended to decrease the critical speed 
slightly. This shows that SFD act only as energy dis- 
sipators not as springs in series with the rotor. Note 
that the resonance peaks are not well defined, appearing 
as eroded versions of what they might have been with 
"rigid" supports. The onset of a noise which was taken 
to be that of cavitation was noted at the peaks where the 
"erosion" begins. Comparing the plots for the various 
clearances, it can be seen that larger clearances yielded 
freer motion of the damper journals and thus "softer" 
dampers . 

Larger clearances yielded more favorable rotor 
responses (for this system) as seen in Figures 5.7 and 
5.8. Considering Figures 5.9,5.10, 5.11, the damper 
orbits tended to "lift off" and centralize with decreas- 
ing clearance. In view of this series of figures, larger 
more centralized orbits with respect to the clearance 
circle provided sti f f er damping and less dissipation to 
the system than did offset orbits. 

Increasing the inlet pressure delayed and lessen- 
ed in intensity the cavitation noise noted earlier and. 
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Figure 5.8 





considering Figures 5.12 through 5.16, also decreased the 
dissipation. Again, examining Figure 5.12, cavitation 
tends to ''erode'' the resonance peak. In view of Figures 
5.17 through 5.20 with increasing pressure the damper 
orbits tended to increase in size and slightly in offset 
until a pressure was reached that centralized the jour- 
nals due to the circumferential inlet effect (see Figure 
5.20 for 30 psi). Holmes suggested that less cavitation 
leads to less load carrying capacity (pg 38) which 
implies larger orbit offsets in agreement with these re- 
sults. So, increasing inlet pressure softens the damper 
but decreases cavitation and energy dissipation. This 
suggests that cavitation is a powerful dissipating device 
since less cavitation yielded less dissipation despite 
the resulting larger orbits and offsets (which tended to 
increase dissipation). 

Figures 5.21 through 5.24 clearly show that rais- 
ing the oil temperature improved the vibrational response 
of the system by increasing dissipation. Increasing the 
temperature lowers the viscosity which allows for freer 
motion (softening) of the dampers. Figures 5.25, 5.26, 
show that with higher oil temperatures the elliptical 
damper orbits are elongated and the offset is increased 
coinciding as before with increased dissipation. 

Figure 5.27 shows damper orbits at three differ- 
ent speeds. As the speed was increased through the 





Midshaft Displacement vs. Rotor Speed 
Groove Feed, 80°F, C = .006", u = 2 gm 


1 P,„ • 2 PSi 



Fiaure 5.12 
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Midshaft Displacement vs. Rotor Speed 
Groove Feed, SO^F, C = .010", u = 3 gm 



Figure 5.14 
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Midshaft Displacement vs. Rotor Speed 
Groove Feed, C = .020°, 80°F, u = 1 gm 



Figure 5.15 


PHASE 
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Midshaft Phase vs. Rotor Speed 

Groove Feed, C ? .020", 80®F, u = 1 .0 gm 

1 Pin = 2 psi 

2 P.^ = 20 psi 



Figure 5.16 



Damper Orbits Relative to Clearance Circle 
Groove Feed, 120®F, u = 3.0 gm, 9000 RPM 



Figure 5.17 


Damper Orbits Relative to Clearance Circle 
Groove Feed, 80®F, u = 3.0 gm, 9000 RPM 



Figure 5.18 
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Midshaft Phase vs. Rotor Speed 

Groove Feed, C = .020", = 2pgi , u = 1.62 gm 

1 

2 120®F 



Figure 5.21 
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Midshaft Displacement vs. Rotor Speed 
Groove Feed, C = .010", = 2 psi, u = 3 gm 



Figure 5.23 





Figure 5.24 
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critical speed, the amplitude grew and the orbit "lifted 
off" of the clearance circle with the higher vibration. 
This figure also portrays the phase shift depicted by the 
reference gap in the orbit tracings. The rotor operated 
clockwise while the gap moved counterclockwise with in- 
creasing speed. 

SFD With End Seals 

Figure 5.28 and 5.29 show the effects on rotor 
response of varying the unbalance with low inlet port 
pressure. The Bode plots are similar to those for SFD 
without end seals operating at low inlet pressure. Cavi- 
tation noise was again noted beginning at the speed where 
the peaks became eroded. As before, increasing unbalance 
created greater excitation and thus more energy dissipa- 
tion. 

With the inlet port at the bottoms of the damper 
"lift off" was aided by the inlet pressure. Figure 5.30 
depicts lift off due to increasing port pressure. 

Increasing the inlet pressure again decreased 
cavitation and lessened energy dissipation (Figures 5.31 
through 5.35) despite creating larger orbits and offsets 
(Figure 5.30). As with SFD without end seals, larger 
clearances or higher temperatures yielded softer damping 
and better performance for this system. 

With low inlet pressure the dampers lifted off 
with speed as with grooved inlets. With higher pressures 









Port Feed, 80°F, = 

1 U = .0166 (.5 gm) 

2 U = .0331 (1 gm) 

3 U = .0993 (3 gm) 
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Figure 5.31 
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Port Feed, C = .02 
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X 1000 RFM 


Figure 5.32 




9 



1 Pin = 15 psi 

2 Pin = 23 psi 

3 Pin = 40 psi 



X ICXXDFFM 


Figure 5.33 


Midshaft Displacement vs. Rotor Speed 









Figure 5.35 









the damper journals lifted off to the top of the clear- 
ance circle at low speed then the critical speed excita- 
tion forced a "lift off" downward as shown in Figure 5.36. 
A medium pressure could cause a lift off with speed that 
centered the orbit at the critical speed as in Figure 
5.37. This is noteworthy since some analysts in the past 
have assumed circular centered orbits near the critical 
speed for ease of analytic solution. 



Damper Orbits Relative to Clearance Circle 


Port Feed, 80®F, u = 3 gm, = 43 psi 



Figure 5.36 
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Damper Orbits Relative to Clearance Circle 
Port Feed, 80®F, u = 2 gm, = 22 psi 



Figure 5.37 




CHAPTER 6 


COMPARISONS WITH A SIMPLIFIED COMPUTER MODEL 

The rotor was modeled using an available code, 
assuming the rotor to be rigid (14). Of course, the 
rotor was not rigid at the speeds which the experimental 
results were obtained. Thus, the damper excitation was 
not necessarily the same and direct comparisons could not 
be made. However, qualitative comparisons concerning the 
general character of the orbits were found to be useful. 
All computed cases presented were run at 10,000 RPM, with 
3 gm unbalance and low inlet pressure using a short bear- 
ing approximation of the Reynolds equation. 

Considering Figures 6.1, 6.2, 6.3, with increas- 
ing viscosity, the orbits' amplitudes and offsets de- 
creased, showing the damper to be growing stiffen. This 
agrees with the experimental results (see Figure 5.25, 
5.26). Also, as the offset decreased, the orbits became 
more circular, which also agress with experimental evi- 
dence. All of the computer produced orbits (clockwise 
rotations) are left of center (e.g. see Figure 6.2) and 
approach the center line (y axis) with decreasing offset. 
Comparing with Figures 5.9, 5.10, 5.11, this is also the 
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Damper Orbits Relative to Clearance Circle 

Groove Feed, Low P._ u = 3.2 gm, 10,000 RPM 

I n j 

1 y = 10"® Reyns. (160®F) 

2 y = 5x10"® Reyns. (100“F) 

3 y = 10x10"® Reyns. (80°F) 



Figure 6.1 
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Damper Orbit Relative to Clearance Circle 
Groove Feed, Low P^.^, u = 3.2 gm, 10,000 RPM 

1 10“® Reyns. (160®F) 

2 10x10"® Reyns. (SO^F) 



Figure 6.3 
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the general rule for the experimentally obtained orbits. 
From this last comparison it can be seen that increasing 
the clearance in both cases produced larger orbits and 
offsets and thus a softer damper. 




CHAPTER 7 


SUMMARY 

In previous analysis a nondi mensi onal damping 
parameter was typically used, originating from the Rey- 
nolds equation, that was proportional to u/C^ (see page 

39), This parameter was used to express the stiffness of 

the damper. As has been shown, experimental evidence 
tends to support this. 

Cookson and Kossa used two nondimensi onal para- 
meters, a bearing parameter (F) and a gravity parameter 
(W), for SFD design and suggested values of these to 
smooth critical speed transition. Their suggestions of 
^ i. • i'^plies that C .004 inches for this system. 

With C = .01 inches, the suggested value of ¥ = .1, re- 

quires the viscosity to be decreased by a factor of five 
from the room temperature value or the temperature to be 
raised to 150®F. This tends to agree with experimental 
evidence. However, for C = .02 inches, an increase in 
viscosity is recommended which is not in agreement. 

In order to eliminate some difficulties that were 
encountered a recommended improvement is the use of 
tapered damper inserts. The present damper inserts were 
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machined to be easily removed from the bearing and 
interchanged, not requiring a pressed fit. This requir- 
ed a radial clearance of one to two mils. This annular 
clearance was filled with oil during damper operation and 
acted as a second SFD in series with the desired damper. 
The idea was that if this secondary SFD had a clearance 
much smaller than the desired damper, it would be much 
stiffer and would not significantly offset the perfor- 
mance of the damper. However, with low inlet pressure 
this secondary damper might have been starved of oil thus 
affecting the damper as a whole. Also, this smaller 
clearance added to the measured damper clearance, so that 
for instance the measured 12 mil clearance was probably 
actually closer to 10 mils. This also provided diffi- 
culty in locating the orbit locations within the clear- 
ance circle, possibly explaining the lateral difference 
in orbit positions on Figures 5,25 and 5.26. These dif- 
ficulties could be overcome easily with the use of taper- 
ed inserts that would be lightly pressed into the damper 
bearings, requiring only slight modifications. 

If increased energy dissipation is desired for 
this rotor system, it has been shown that the dampers 
need to be softened. A lighter lubricant is recommended 
to avoid heating the oil or using larger clearances. 

Large clearances are both unusual in practice and 
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difficult to model due to the necessity of involving 
inertial effects ( 16 ). 

The parametric study conducted here, hopefully 
provided some insights into the workings of an SFO, 
especially concerning the effects of cavitation. More 
importantly, the test rig has been shown to be capable of 
providing useful data so that with the above mentioned 
improvement, data can be taken with a more specific pur- 
pose to be compared with computer analysis. 
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APPENDIX A 


DESIGN CALCULATIONS 

A rotor shaft was to be designed to operated on 
the existing Bent ly-Nevada test rig platform. This con- 
strained the bearing to bearing shaft length to 18" or 
less. The use of proximity probes to measure displace- 
ment required the material to be carbon steel. Bearing 
and damper considerations suggested a diameter of 13/16". 
With the rotor disks positioned in the center of the 
shaft, it was desired that the rotor would have its first 
bending critical speed near the end of the operating 
speed range (0 to 12,000 RPM), The critical speed was 
chosen to be about 9000 RPM. Given the above considera- 
tions, the bearing to bearing length of the shaft was to 
be determined. 

The stiffness (k) at the center of a simply sup- 
ported beam is 

48EI 

< = -T3- 

where E is the modulus of elasticity and I is the section 
modulus. An estimate of the first natural frequency of 
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1/2 

the shaft is then o»c = /<7ni where 
at the center of a massless shaft, 
equations and solving for L, 

f48En i/3 


L = 


mo)' 


m is a lumped mass 
Combining the above 


Using the mass of the two rotor disks (2.5 lb) as the 
lumped mass and the above mentioned critical speed and 
dimensions, this equation yields 17.5". Calculating the 
mass of a shaft of this length (3 1b) and lumping half of 
this mass in the center in addition to the disks, the 
above equation then yields 15.0". This was taken as the 
final estimate since the design allowed for some adjust- 
ment room. This estimate proved to be satisfactory, 
yielding a critical speed of about 9200 RPM. 

Using the above equations along with the experi- 
mentally measured critical speed, the above method was 
reversed yielding the bearing and disk masses necessary 
to simulate the shaft as massless for operation through 
the first critical speed. These values were helpful for 
data analysis. Solving for the disk mass. 


m - 48EI 
'"0 - 


This equation yields mg = 3,992 1b using the appropriate 
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design parameters. The total rotor mass was measured as 
6.591 lb. Half the difference of the disk and total 
rotor masses equaled the mass to be lumped at each bear- 
ing station, mg = 1.300 lb. So that the mass ratio (M) 
was 1/3. 

A quill shaft was to be designed to be used as a 
flexible coupling between the rotor and rotor shaft. The 
quill had to be strong enough to transmit the maximum 
torque of the motor, yet flexible in bending so as to 
constrain minimally the translation of the rotor. Also, 
the first critical speed of the quill was required to be 
well above the operating speed of the rotor, so that the 
dynamic behavior of the quill would not offset that of 
the rotor. 

Considering the dimensions of the rotor, a dia- 
meter of 1/8" seemed reasonable for the quill. The motor 
was rated at 1/10 hp. at 10,000 RPM or a maximum torque 
(T) of .1728 in. -lb. From the equation 

Tr 

^ = -J- 

the maximum shear stress was calculated to be 450 psi. 
This was deemed suitable for a steel shaft. 

Given the length of the rotor and the length con- 
straint of the test rig base, a maximum quill length of 
6" was determined. The first bending natural frequency 



ni 


of a clamped beam is given by 
f 


_ (4.731^ fElV/2 

2nL^/^ [m J 


For a 6" steel shaft with a 1/8" diameter, this equation 
yields f = 37,600 RPM which is over three times the 
operating speed of the rotor. This design was then con- 
sidered satisfactory. 
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Figure B-5. Wiring Schematic 
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Naval Air Propulsion Test Center 
Aeronautical Engine Department 
Trenton, NJ 08628 
Attn: Mr. Robert De Lucia 
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12 New England Executive Park 
Burlington, MA 01803 
Attn: Mr. Robert Berman 
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Office of Aviation Safety, FOB lOA 
800 Independence Ave., S.W. 

Washington, DC 20591 

Attn: Mr. John H. Enders 1 

FAA, ARD-520 

2100 2nd Street, S.W. 

Washington, DC 20591 

Attn: Comnander John J. Shea 1 

National Transportation Safety Board 
800 Independence Ave., S.W. 

Washington, DC 20594 

Attn: Mr. Edward P. Wizniak TE-20 1 
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and Engineering Science 
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Attn: R. Loewy 1 

Cleveland State University 
Department of Civil Engineering 
Cleveland, OH 44115 
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P.O. Box 3999 
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Attn: Dr. R. J. Bristow 
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P.O. Box 3707 

Seattle, WA 98124 

Attn: Dr. Ralph B. McCormick 
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P.O. Box 3707 

Seattle, WA 98124 

Attn: Mr. David T. Powell 

Boeing Commercial Airplane Company 

P.O. Box 3707 

Seattle, WA 98124 

Attn: Dr. John H. Gerstle 

Boeing Company 
Wichita, KA 
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Lambert Field, MO 63166 
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P.O. Box 1128 
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Attn: Dr. L. Beitch 

Dr. M. Roberts 
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Technology Center 
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Hawthorne, CA 90250 
Attn : Library 

North American Rockwell, Inc. 
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Downey, CA 90241 
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Worcester, MA 01606 
Attn: Mr. George E. Buron 

Norton Company 
1 New Bond Street 
Industrial Ceramics Division 
Worcester, MA 01606 
Attn: Mr. Paul B. Gardner 

United Aircraft Corporation 
Pratt 6 IVhitney Group 
Government Products Division 
P.O. Box B2691 
West Palm Beach, FL 33402 
Attn : Library 

R . A . Marmo 1 

United Aircraft Corporation 
Pratt § Whitney Aircraft Group 
400 Main Street 
East Hartford, CT 06108 
Attn: Library 

R. Liss 
D. H. Hibner 
C. Platt 
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Attn: F. Sisto 

A. T. Chang 
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Latham, NY 

Attn: M. S. Darlow 
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Northway 10, Executive Park 
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P.O. Box 1103 
Huntsville, AL 18908 
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Attn: R. H. MacNeal 

MARC Analysis Research Corporation 
260 Sheridan Avenue, Suite 314 
Palo Alto, CA 94306 
Attn: P. V. Marcel 

United Technologies Research Center 
East Hartford, CT 06108 
Attn: Dr. A. Dennis 

Georgia Institute of Technology 
School of Civil Engineering 
Atlanta, GA 30332 
Attn: S. N. Atluri 

Georgia Institute of Technology 
225 North Avenue 
Atlanta, GA 30332 
Attn: G. J. Simitsis 

Lavrrence Livermore Laboratory 
P.O. Box 808, L-421 
Livermore, CA 94550 
Attn: M. L. Wilkins 

Lehigh University Institute of Fracture 
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Bethlehem, PA 18015 
Attn: G. T. McAllister 

Materials Science Corporation 
1777 Walton Road 
Blue Bell, PA 19422 
Attn: W. B. Rosen 

National Bureau of Standards 
Engineering Mechanics Section 
Washington, DC 20234 
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Purdue University 
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West Lafayetee, IN 47907 
Attn: C. T. Sun 

University of Dayton Research Institute 
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Attn: F. K. Bogner 
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Attn: R. H. Gallagher 

H. Kamel 

University of California 
Department of Civil Engineering 
Berkeley, CA 94720 
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